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We propose and investigate bounds on quantum process fidelity of quantum filters, i.e. proba¬ 
bilistic quantum operations represented by a single Kraus operator K. These bounds generalize 
the Hofmann bounds on quantum process fidelity of unitary operations [H.F. Hofmann, Phys. Rev. 
Lett. 94, 160504 (2005)], and are based on probing the quantum filter by pure states forming two 
mutually unbiased bases. Determination of these bounds therefore requires much less measurements 
than full quantum process tomography. We find that it is particularly suitable to construct one of 
the probe basis from the right eigenstates of K, because in this case the bounds are tight in the 
sense that if the actual filter coincides with the ideal one then both the lower and upper bounds 
are equal to one. We theoretically investigate application of these bounds to a two-qubit optical 
quantum filter formed by interference of two photons on a partially polarizing beam splitter. For 
experimentally convenient choice of factorized input states and measurements we study the tight¬ 
ness of the bounds. We show that more stringent bounds can be obtained by more sophisticated 
processing of the data using convex optimization and we compare our methods for different choice 
of the input probe states. 

PACS numbers: 03.67.-a, 03.65.Wj 


I. INTRODUCTION 

Characterization of quatum processes represents an in¬ 
dispensable tool for testing, optimization, and bench¬ 
marking of quantum information processing devices. A 
complete characterization of a quantum device can be 
provided by quantum process tomography mm, but this 
technique becomes very time consuming with increasing 
complexity of the device, unless special procedure such 
as compressed sensing can be applied HE], because the 
number of parameters that need to be estimated scales 
exponentially with the number of qubits. 

Instead of full quantum tomography we may just at¬ 
tempt to obtain an indication how close we are to the 
target operation, as quantified by the quantum process 
fidelity. Monte Carlo sampling has been proposed for 
efficient estimation of fidelity of multiqubit states and 
operations with resources scaling polynomially with the 
estimation precision 0 H]. Yet another experimentally 
appealing option is represented by the Hofmann bounds 
on quantum process fidelity [7j. In this approach, a 
lower and upper bound on fidelity with a unitary op¬ 
eration is determined from measurements of average out¬ 
put state fidelities for input states forming two mutually 
unbiased bases. This latter approach is particularly ef¬ 
ficient for characterization of few-qubit operations and 
it can be used e.g. for preliminary benchmarking of a 
quantum device during its design and optimization be¬ 
fore a more complete characterization is carried out at 
the optimal operating point. During recent years, the 
Hofmann bound has been successfully utilized for exper¬ 
imental characterization of various two-qubit and three- 
qubit quantum gates [SHIS]. 

The Hofmann bound was designated to provide bounds 
on quantum process fidelity with a deterministic unitary 
operation. Here we generalize this technique and propose 


and investigate Hofmann-like quantum process fidelity 
bounds for special kind of probabilistic quantum opera¬ 
tions called quantum filters. Quantum filters are com¬ 
pletely positive trace decreasing maps that can be rep¬ 
resented by a single Kraus operator K. Quantum filters 
form an important tool in many branches of quantum in¬ 
formation science and beyond and they find applications 
e.g. in quantum state engineering, entanglement distilla¬ 
tion |I6H18] . or quantum state discrimination [SI HOI- 

Our derivation of the generalized Hofmann bounds for 
quantum filters is based on operator inequalities that are 
at the heart of the original Hofmann bound. In contrast 
to unitary operations, where measurement of state fideli¬ 
ties for two complementary bases is sufficient, in case 
of quantum filters one generally needs to perform an 
additional set of measurements, which essentially char¬ 
acterizes the performance of the filter in a basis of its 
eigenstates. The number of measurements can be kept 
the same as for unitary operations provided that one of 
the two input bases is formed by the right eigenstates of 
K. In this case it can be also proved that the result¬ 
ing lower and upper bounds on quantum process fidelity 
are tight in the sense that for a perfect filter the bounds 
are always equal to 1. The probabilistic and non-unitary 
nature of the quantum filters thus leads to a symmetry 
breaking and occurrence of a preferred set of probe input 
states. We explicitly consider two bases connected via 
Fourier transform and also two n-qubit bases connected 
by Hadamard transform on each qubit. 

As an illustration, we theoretically investigate charac¬ 
terization of a two-qubit optical quantum filter formed 
by interference of two photons on a partially polarizing 
beam splitter followed by post-selection of detection of 
a single photon at each output port of the beam split¬ 
ter. This filter is utilized in various linear optical quan¬ 
tum information processing devices such as linear optical 
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quantum gates [i m mHia. We consider experimen¬ 
tally convenient choice of input probe states for which 
the required output state fidelities can be directly deter¬ 
mined by product single-qubit measurements. We show 
that as a consequence of this basis choice the resulting 
upper and lower bounds are not tight. We numerically 
find the ultimate upper and lower bounds for the same 
data using the semidefinite programming approach. In 
this way we illustrate that for the considered quantum 
filter and the available data more stringent bounds can 
be obtained by more sophisticated processing of the data 
and we compare the two methods also for different choice 
of basis states. 

The rest of the paper is organized as follows. In section 
[IT] we review the original Hofmann bound on fidelity of 
a deterministic process with a fixed unitary process and 
we generalize this bound to quantum filters. In section 
|III| we study quantum filter formed by interference of two 
photons on a partially polarizing beam splitter and sec¬ 
tion IV contains the conclusions. Finally, the Appendix 
contains a proof of an alternative lower bound for n-qubit 
filters and some technical derivations. 


II. FIDELITY BOUNDS FOR QUANTUM 
FILTERS 

In general any quantum operation can be repre¬ 
sented using Choi-Jamiolkowski isomorphism [24l[23 by 
positive-semidefinite operators. In this way a quantum 
filter T ■. p ^ KpK^ is represented by an operator 
XT = I, where \ujk) = I ® K\uj), 

l‘") = 4Ele.)|e.) (1) 

V « 

denotes a maximally entangled state and is an 

orthonormal basis of a d-dimensional Hilbert space Hd- 
Fidelity between the actually implemented quantum op¬ 
eration X s-nd the ideal quantum filter xt can be defined 
via state fidelity between normalized Choi operators [50] j 

p ^ Tr(x Xt) .21 

Tr(x)Tr(x^)- ^ ^ 

Our aim is to propose a procedure that would lower and 
upper bound this fidelity based on measured state fideli¬ 
ties of output states with respect to the ideal output for 
several input states. To do this we first review the orig¬ 
inal Hofmann bound for deterministic operations in d 
dimensions |7| and then we generalize it so that it will 
become applicable to non-unitary quantum filters. 

Suppose {\fk)}t=i are two orthonormal 

bases of Hd that are mutually related by discrete Fourier 
transform, i.e. 


\fk) 






(3) 


We denote by |e°“‘) = U\ej), |/™*) = U\fk) the ideal 
output states of a unitary transformation U. Density 
matrices pj and of output states produced by the ac¬ 
tually implemented operation x from the input states 
|ej), l/fc) can be expressed as 

Pj = dTiinix |ej)(ejf 0 7), 

a = dTri„(x|/fc)(/fcr0/), (4) 

where the transposition is taken with respect to the basis 
|ej) used in the Choi-Jamiolkowski isomorphism. The 
average output state fidelities for the two sets of probe 
states are defined as follows 

'"I = 

^2 = ^E(/r‘iai/r)- (5) 

The Hofman lower bound on quantum process fidelity |7] 

F>Fi+F2-1 (6) 


can be proved as follows. For deterministic transforma¬ 
tions Tr(x) = 1, hence the bound (§ is equivalent to 

d 

Tr(x \u;u){u^u\) > E |e,)(e,r ® |er)(cri) 

i=i 

+ETr(xi/fc)(/.r®i/r‘)(/r‘i) 

k^l 

-Tr(x/0/). (7) 

The validity of the inequality 0 would be guaranteed 
by showing the positive-semidefiniteness of an operator 
X = I®URI®U\ where 

d 

i?=|w)(a;| -^|ej)(ejf 0 |ej)(ej| 
i=i 

d 

( 8 ) 

k=l 

If X > 0 then Tr(xX) > 0 due to positivity of x > 
0, which implies the inequality Q. As we show in the 
appendix [A| the operator R can be rewritten (in the term 
to term fashion) as 


d d 

R =|wii)(wii| - ^ ^ |wifc)(a;ifc| 

j=l k=l 

d 

+ 'y ^ \^jk) {^jk\: (9) 

j,k=l 

where {|wjfe)} is an orthonormal basis of maximally en¬ 
tangled states in d dimensions. 
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From the above equation it is clear that i? > 0, which 
implies X > 0 and this proves the original Hofmann 
bound. In a similar fashion, the upper bounds on quan¬ 
tum process fidelity F < Fi, F < F 2 can be derived from 
the following two operator inequalities 

d d 

|wii)(a;ii| < |a;ifc)(a;ifc|. 

j=l k=l 

( 10 ) 

Next, we will use the operator R to derive lower and 
upper bound on the fidelity of quantum filters. Let us 
multiply Rhj I ® K from the left and by / 0 from 
the right. Using Eq. ^ we obtain 

d 

\^k){ujk\ - \ej){eY ® K\ej){ej\K'^ 
i=i 

d 

k=l 

where the inequality follows from R > 0. By taking the 
trace with x Ed- 0 can be rewritten as 

1 ^ ^ 

Tr(x \u}k){u}k\) - ^^Tr(iF|ej)(ej|iF^ pj) 
i=i 

1 

- ik) + TriKK^ H) > 0, 

( 12 ) 

where pj and defined in Eq. Q are the unnormalized 
output states of a probabilistic operation y corresponding 
to pure input states jcj) and \fk), respectively, and fl = 
dTrin(x ® I) is the unnormalized output state for 

a maximally mixed input state. To obtain a lower bound 
on fidelity of a quantum filter K, we divide the inequality 
(121 by Tr(x) Tr(xjF) and rewrite the resulting expression 
such that it contains normalized overlaps, 

d d 

E >^Pj(ej|ft|ej) +Y^k{fk\ik\fk) - A Tr^KK'^ Cl). 

j=l k=l 

(13) 


Here pj = pj/Tv{pj), Ik = Cfc/Tr(^fe), and Cl = 
r2/Tr(H), are normalized output states of x ^md 

K\fk) 

VieAK^KleY V{h\K^K\fk)^ 

denote the normalized output states of the ideal filter K. 
The weight factors read 

p, = APj{ej\K^K\e,), g,. = AQk{fk\K^ K\fk), 

(15) 

where A = d/TT{K'^K) = 1/Tr(xjr), and Pj = 

and Qk = dTrix) the relative success proba¬ 

bilities of operation x for input basis states \ej) and 


l/fc), respectively. These relative probabilities satisfy 
Pj — 12k Qk = 1- 

Formula (131 generalizes the Hofmann lower bound ([^ 
to quantum filters and represents one of our main re¬ 
sults. We can see that the hdelity of a quantum filter 
is lower bounded by an expression which contains two 
weighted sums of the output state fidelities {ej\pj\ej) and 
{fk\^k\fk), which generalizes the average state fidelities 
El and F 2 appearing in the original Hofmann bound. The 
last term on the right-hand side of inequality (131 pro¬ 
vides a generalization of the factor — 1 to quantum filters 
and depends both on the ideal filter K and also on the 
actual operation x through Cl. It follows from Eq. (12) 
that the summation in Eq. should be performed 

only over those terms for which the overlap {ej\K^K\ej) 
or {fk\K^K\fk) is nonzero. This also ensures that the 
normalized output states (14) are well defined. 

A most straightforward way to experimentally deter¬ 
mine the output state fidelities and relative success prob¬ 
abilities Pj and Qk consists in measuring the output state 
Pj {^k) in a basis including the corresponding output 
state \ej) {\fk)) produced by an ideal filter. For quantum 
filters {\ej)}'j^i and {\ fk)}k=i generally do not form a ba¬ 
sis which means that a separate measurement basis must 
be set for each probe state. Besides testing the unknown 
quantum transformation with 2d input states {|ej)}^^i, 
{|/fc)}fc=ij ■^e also need to determine the term Fv{KK^Cl) 
by some measurements. To construct a suitable measure¬ 
ment, consider a singular value decomposition of AT, 


d 

1=1 

where the left and right eigenvectors ju;), \wi) form two 
orthonormal bases and the non-negative singular values 
were chosen in the form -y/A to simplify subsequent for¬ 
mulas. As a consequence the positive-semidefinite oper¬ 
ators K, KK^ have the following spectral decomposi¬ 
tions 


d d 

K^K = YMwi){wiI KK^ = YMvi){vi\. (17) 
1=1 1=1 

In principle, we can determine Tr{KK^Cl) from suit¬ 
able measurements on any d input states forming an 
orthonormal basis (e.g. vectors \uj)). Let Q = 
(iTrin(x \uj){uj\'^ 0 I) and (j = Ci/Tr(Cj) denote the 
unnormalized and normalized output state correspond¬ 
ing to the input state \uj), and similarly as before we 
define the relative success probability for this input state 
as Rj = The term Tr{KK^ H) can then be ex¬ 

pressed as 


Tr(EEt Cl)=Y XiRj{vi\Q\vi). (18) 

3,1 = 1 
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The relative success probabilities Rj as well as the over¬ 
laps of output states C,j with |u;) can be determined by 
measuring each output state C,j in the basis formed by 
the eigenstates \vi). 

At this point it is useful to realize that the experimen¬ 
tal effort can be kept the same as for K being unitary 
at the price of a suitable choice of basis \ej) and con¬ 
sequently l/fc). Thus, if we choose \ej) = \wj) and also 
l^i) = \wj) then 


,) = 


K\wj) 


\/{wj\K^ K\wj) 


= \vj) 


(19) 


and the data for input states \ej) measured after the filter 
in basis je^) can be used to determine the last term of 
Eq. (13), i.e. Ci = Pj^ Rj = Pj Vj- After some algebra 
we fiiMAhat the lower bound now equals to 


a a X 

F - E t (1 - W{ei\p,\~ei), 

k^l j,l^l ^ 


( 20 ) 


where A = (X)j=i used the identity 

{fk\K^K\fk) = A, which holds since the two bases \ej) 
and l/fc) are related by quantum Fourier transform. 

An important property of the lower bound is its tight¬ 
ness. Especially, if the implemented transformation is 
the desired one, then fidelity F = \ and we want our 
lower bound to attain the value 1 as well. If the im¬ 
plementation of the filter is perfect, then (/fc|^fc|/fc) = 1 
and {ei\pj\ei) = 5jk- If we insert these expressions into 
Eq. ( |20| ) then we get F > 1 which confirms that the lower 
bound (20) is tight for any quantum filter K. Note that 


this tightness is achieved due to the special choice of the 
probe states, where |ej) coincide with the right eigenvec¬ 
tors of K. For other choices of the probe states the lower 
bound (13) is generally not tight and can be strictly lower 


than 1 even for a perfect filter. This should be contrasted 
with the original Hofmann bound which always at¬ 
tains value 1 if the target unitary U is implemented per¬ 
fectly, irrespective of the choice of the two probe bases. 
The nonunitarity of the filter K thus leads to a symme¬ 
try breaking and emergence of preferred states suitable 
for benchmarking of the filter. 

In a very similar fashion as above also a pair of upper 
bounds can be derived 


F <^Pj{ej\pj\ej), F < ^9fc(/fc|/fe|/fc), (21) 


i=i 


k=l 


that are equivalent to the inequalities (10). For the spe¬ 
cial choice of probe states \ej) = \wj) the upper bounds 
simplify to 




i\pj\vj), F < ^Qfc(/fc|/fc|/fc). 


/c=l 


(23) 


Similarly as the lower bound (20), also the upper bounds 


(23) are tight in the sense that they yield F < 1 if the 
filter is implemented perfectly. 

From the application point of view, the n-qubit sys¬ 
tems whose Hilbert space is endowed with tensor product 
structure and d = 2" are particularly relevant. In this 
case a natural choice of \ej) could be the computational 
basis formed by tensor products of single-qubit states |0) 
and |1). By its construction, discrete quantum Fourier 
transform of the n-qubit computational basis states leads 
to product n-qubit states |/fc). This is good for experi¬ 
ments, since preparation of product states is often much 
simpler than preparation of entangled states. Unfortu¬ 
nately, for most quantum filters at least some (ideal) 
output quantum states for the above inputs are entan¬ 
gled. For that reason it might be useful to study also 
other pairs of input bases, which could be tested more 
easily. One such combination can be the computational 
basis and the Hadamard basis formed by tensor products 
of states |±) = (|0) ± \l))/^/2. As we show in the Ap- 
pendix[^the above derived lower and upper bounds hold 
also for this latter setting. 


III. TWO-QUBIT QUANTUM FILTERS 


The Hofmann bound (|^ proved particularly useful for 
characterization of various linear optical quantum gates 
[HHumiin]- As a case study, we therefore investigate 
here a characterization of a quantum filter acting on po¬ 
larization state of two photons. The filtering is realized 
by interference of the photons on a partially polarizing 
beam splitter (PPBS) followed by postselection of events 
when a single photon is observed at each output port of 
the beam splitter. In the basis of vertical and horizon¬ 
tal polarizations, the resulting two-qubit quantum filter 
reads 


K = 



0 

tutv 

-rnrv 

0 


0 

-thTy 

tntv 

0 



(24) 


if we start from inequalities 

d 

|w)(w| <^\ej){ej\^ ®\ej){ej\, 

d 

M(a;|<^|/fc)(/fcf ®|/fc)(/fc|, (22) 

fc=i 


where t^r, ty and rn, ry denote the amplitude transmit- 
tances and reflectances of PPBS for horizontal and verti¬ 
cal polarizations, respectively. We assume that the trans- 
mittances and reflectances are real and that the beam 
splitter is lossless, hence t’j + Vj = 1. 

We will investigate further only the case t// = 1, 
since this element is often used in optical quantum in¬ 
formation processing experiments and its fidelity should 
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FIG. 1. A lower bound (131 on fidelity F of a two-qubit 
quantum filter K is plotted in dependence on the intensity 
transmittance Tv = ty of the partially polarizing beam split¬ 
ter. Ideal implementation of the filter is assumed, hence the 
true fidelity F = 1 and is depicted by the blue dashed line. 
The gap between the two lines illustrates the tightness of the 
bound. 






FIG. 2. The red solid lines show the dependence of the lower 
bound and fidelity of a quantum filter K on the actual in¬ 
tensity transmittance Tv = ty of the PPBS for four different 
target transmittances 0.1 (a), 0.25 (b), 0.5 (c), and 0.75 (d). 
The true fidelity of the filter is depicted by blue dashed lines. 


where 


be assessed. In this case the filter becomes diagonal, 
K = A\ag{l,tv,tv,2ty - 1}. Choosing as 

the computational basis in which the filter is diagonal 
would for probe states \fk) lead to measurements on out¬ 
put states in an entangled basis, which is problematic in 
many optical experimental setups. Instead, we will in¬ 
troduce alternative probe states that require just prepa¬ 
rations and measurements in product bases. The idea is 
to employ the following pair of bases. 


|ei) = |0)|+), 

|e2) = |0)|-), 

|e3) = |l)|+), 
|e4) = |l)|-), 


l/i) = l+)|0), 

1 / 2 ) = l+)|l), 

l/3) = |-)|0), 
l/4) = |-)|l). 


(25) 


The choice of these probe states is motivated by previ¬ 
ous experiments, where bounds on fidelity of a quantum 
controlled-NOT gate and controlled-Z gate were deter¬ 
mined [HHIIllIi]- The two bases are related via a 
Hadamard transform on each qubit, \fj) = H 0 H\ej), 
and this relation together with the factorized form of the 
basis states \ej) and \fk) ensures that the lower and up¬ 
per fidelity bounds (13) and (21) are applicable, c.f. also 


Appendix B. The practical advantage of the probe states 
(25) is that the filter K maps them on product states. 


|ei) = |0)|a+), 
lea) = |0)|a_), 
|e3) = |l)IM, 

|e4) = |l)|&-), 


1/1) = |a-H)| 0 ), 

1 / 2 ) = |6+)|1), 

l/3) = |a-)|0), 

1 / 4 ) = I&-)|1), 


|a±) = 


\b±) = 


| 0 )± 


tv 


\/l + b-v \/l + b 


^|1), 


W 


^ t ' i , + (24 - 1)2 


| 0 )± 


2ty 1 


■\/tv + (2ty — 1)2 


|1)- 


Since the filter K is not diagonal in any of the two 
probe bases (25), additional measurements are required 
to estimate the term Tr{KK^(l). A natural choice is 
to employ the computational basis states as additional 
probes \uj) and measure the resulting output states in 
the computational basis. Mathematically, we have \uj) = 
\wj) = \vj) and 

|ui) = |00), \V2) = |01), \V3) = |10), |z;4) = 111). 

(26) 

The above outlined procedure to test the unknown quan¬ 
tum filter requires 3d = 12 product input states and each 
output state needs to be measured in a single product 
basis. For comparison, full quantum process tomogra¬ 
phy of a two-qubit quantum filter m would typically 
involve about 144 different combinations of input states 
and output measurements. Thus, characterization of the 
filter via the generalized Hofmann bounds requires much 
less resources. 

Let us now discuss the tightness of the lower bound 


(13) for our choice of the filter and the probe states. 


A detailed analysis reveals that in the present case this 
bound is tight for the ideal quantum filter, i.e. equal to 
one, only if Ty = ^ or Ty = 1, where Ty = ty. Other¬ 
wise there is a gap, which is smaller than 1% if the desired 
PPBS has Tv > \, see Fig . |h We also considered the 
situation where the actual value of the transmittance Ty 
differs from the desired one. The dependence of the lower 
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FIG. 3. Graphs of upper and lower bounds on the fidelity of two-qubit quantum filters. Upper and lower most dots correspond 
to analytical bounds for randomly generated filters constructed as a mixture of a fixed ideal filter with transmittance tv specified 
in each panel and a randomly chosen filter. The dots just above/below the solid line correspond to the bounds obtained from 
the same data, but using semidefinite programming approach. The results in the upper row correspond to probing with states 
(25 I for which measurements in product basis are sufficient. For comparison, the three graphs in the lower row were obtained 
for a different set of probe states with \ej) chosen to be the computational basis where the ideal filter K is diagonal. This latter 
choice of probe states requires that some of the output states are measured in entangled basis. 


bound on the actual value of Ty is plotted in Fig. [^for 
four different desired values of Ty. The results indicate 
that the proposed lower bound works well for Ty > f. 
However, for Ty < ^ the bound quickly becomes too 
loose to be useful. 

This opens the question what are the best upper and 
lower bounds on fidelity given certain data from exper¬ 
iment. Such a question can be precisely answered by 
semidefinite programming |28) . since the fidelity and the 
constraints given by the measured data are linear in 
Thus, we might write 

Fub = max A Tr(H|a;iy) |) 

A>0,Tr(T) = l 
\/kTT{AMk)=rk 

Flb= min A Tr(H|wiy)(a;iy|) (27) 

A>0, Tr(A) = l 
'ikTv(AMk)=rk 


where the matrices Mj, and parameters capture the 
linear constraints provided by the data. To make our 
formulation of these constraints sufficiently general, con¬ 
sider a set of input probe states forming a 

basis, where each output state is measured in a generally 
different basis {\njk)}k=i- fjk denote the measured 
frequencies which sample the theoretical probabilities 

Pjk = dTv{x \mj){mj\^ ® \njk){njk\)- (28) 

Using the identity k=iPjk = '^Tr(x) we can express 
the set of constraints on A imposed by the data fjk as 


follows. 


Tr(H \mj){mj\^ ® \njk){njk\) = 


fjk 


Si,m=l fl^ 


(29) 


In our present case, each probe basis |ej), \fk), and \vi) 
provides a set of 15 linearly independent constraints. 

To test the performance of this approach, we used it 
to find the best upper and lower bounds from the data 
that we generated for 1000 randomly chosen quantum 
operations y that were constructed as a random mixture 
of an ideal quantum filter K and some other randomly 
chosen filter K' . The optimization (27) was performed 


numerically with the use of CVX, a package for specify¬ 
ing and solving convex programs [HEn]. For a given 
random quantum filter y we also compute the analytical 
lower and upper fidelity bounds according to our proce¬ 
dure and the true fidelity. In figure|^we plot the obtained 
results as a function of the true fidelity, hence the graph 
of the true hdelity forms a straight line with a unit slope. 
For comparison, we consider probing with two sets of in¬ 
put states. The first set is specified by Eq. (25) and re¬ 


quires only product measurements on output states, but 
the analytical bounds are not tight. The second set of 
probe states is constructed such that it leads to a tight 
analytical bounds (20) and (23) at the price of the re¬ 


quirement of measurements in entangled basis for some 
of the output states. Interestingly, the approach based 
on semidefinite programming provides in both cases tight 
bounds, i.e. for ideal implementation of the filter the up- 
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per and lower bounds coincide and J^lb = -Fub = 1- 


IV. CONCLUSIONS 

In summary, we designed and analyzed bounds on 
quantum process fidelity F of a specific type of nonde- 
terministic operations called quantum filters. These op¬ 
erations are mathematically characterized as completely 
positive maps, which can be expressed by a single Kraus 
operator K. Operationally they correspond to opera¬ 
tions which succeed only with a limited probability that 
depends on the input state and they map any pure state 
again to pure state. The proposed bounds represent a 
generalization of the original Hofmann bounds on hdelity 
of unitary transformations [7]. For quantum hlters, the 
average state hdelities are replaced with specific weighted 
averages, and the lower bound contains also an additional 
term that depends both on the desired and the actual op¬ 
eration. As a consequence, in addition to determination 
of relative success rates and output state fidelities for two 
sets of input basis states, further measurements are gen¬ 
erally needed. Nevertheless, we show that the number 
of input states and measurements can be kept the same 
as for unitary operations if one of the two input bases is 
formed by the right eigenstates of K. An important prop¬ 
erty of any bound is its tightness. In particular, for quan¬ 
tum process fidelity bounds we would like to have both 
the upper and lower bounds equal to one if the actual and 
the desired quantum transformation coincide, because in 
that case the fidelity is F = 1. We demonstrated that 
our bound is tight if one set of probe states is formed by 
right eigenstates of K and the other by their quantum 
Fourier transform or by their Hadamard transform. The 
proposed bounds extend the toolbox of efficient methods 
of characterization of quantum operations [Ha ISIH33] 
and provide a method for quick checking of quality of 
quantum filters before their more comprehensive charac¬ 
terization, e.g. by quantum process tomography. 

As an illustration, we have theoretically investigated 
application of the proposed fidelity bounds to character¬ 
ization of a specific two-qubit linear-optical quantum fil¬ 
ter. This hlter is implemented by two-photon interference 
on a partially polarizing beam splitter followed by condi¬ 
tioning on emergence of a single photon at each output 
port of the beam splitter, and it is often utilized in op¬ 
tical quantum information processing with polarization 
encoded qubits. We consider experimentally convenient 
choice of product input probe states for which the re¬ 
quired output state fidelities can be directly determined 
by product single-qubit measurements. It turns out that 
the price to pay for this experimental convenience is that 
the resulting bounds are generally not tight. We com¬ 
pare our analytical bounds with ultimate lower and up¬ 
per bounds that can be obtained from a given experi¬ 
mental data with the help of convex optimization. The 
experimental data represent a set of linear constraints 
and we numerically solve a so-called semidefinite pro¬ 


gram that among all quantum operations satisfying given 
linear constraints finds an operation with minimum and 
maximum overlap with the target quantum filter K. We 
observe that the ultimate lower and upper fidelity bounds 
obtained in this way are tight, i.e. equal to one for perfect 
filters, even if the analytical ones are not. In this way we 
illustrate that for the considered quantum filter and the 
available data more stringent bounds can be obtained by 
more sophisticated data processing. 
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Appendix A: Positivity of operator R 

Let us define an orthonormal basis {\oJjk)}'j u=i max¬ 
imally entangled states via the action of a pair of opera¬ 
tors 

3 3 

on the state |a;) as 

\uj^k) = (A2) 

By definition Iwn) = jw) and 

d ^ d 

E |wij)(wijj = ^ E '’|efc)|efc)(ei|(ez| 

i=l 

d 

<Z\ej}{ej\, (A3) 

because here the transposition is taken with respect to 
the basis \ej). Similarly, 

d ^ d 

E ^ E \'^kej)\ek){ei^j\{ei\. (A4) 

i=i 

On the other hand, 

d 

'^\fk){fk\^ Z) \fk){fk\ = 

k^l 

= 2 E |ej)|ez)(e,„|(e„| 

j,k,l,m,n—l 

= E |ej)lo)(em|(e„|, (A5) 

j,k,l,m,n 
j — l—m—n 




which is clearly equivalent to (A4). Thus the operator R 


defined by Eq. (1^ can be expressed as 


d d 

R =|wii)(a;ii| - ^ ^ \uJik){uJik\ 

i=i 


fc=i 




(A6) 


where we used the identity I ® I = Yl'j k-i \^ 3 k){<^jk\- 
Since we managed to rewrite i? as a sum of projectors, 
R = k =2 \^jk){‘^jk\, this proves that R> 0. 


Appendix B: Alternative lower bound for n-qubit 
filters 


Our goal is to prove positivity of operator R defined 
in Eq. ® for a different pair of orthonormal bases 
{|ej)}^=i, {|/fc)}fc=i- This would allow us to exactly re¬ 
peat the same steps as in the main text and thus we 
could use all the derived lower and upper bounds, but 
for a different choice of \ej), \fk)- Specifically, we shall 
consider systems of n qubits, hence d = 2". Let {\ej)}1j^i 
be the computational basis, where \ej) = \ji)\j 2 ) ■ ■ ■ \jn) 
and jm is the m-th digit in the binary representation of 
number j — 1. Similarly, let {\fk)}k=i be the compu¬ 
tational basis transformed by the Hadamard transform 
(acting as H\jm) = (|0)-|-(—l)'^"*|l))/-\/2) on every qubit, 
\fk) = H\ki) (g) H\k 2 )... (g) H\kn)- Operator R is acting 
on 2n qubits, which are ordered as n qubits related to the 
input state tensored with another n qubits related to the 
output state of the quantum filter. It is useful to divide 
the 2n-qubit Hilbert space on which operator R acts into 
two-qubit subsystems formed by the m-th qubit of the 
input and the m-th qubit of the output. We introduce 
a unitary operator W which groups together the m-th 


input and output qubits [IT] . 

bE|jl . . . 5 Jn) 1^15 ■ ■ • ) ^n) — |jl 5 ^l) ■ ■ ■ \jn} ^n) • (bH) 

In this way the maximally entangled state can be written 
as W\uj) = |$+)i • • • !$+)„, where = ^(|00)±|II)) 
are the Bell states and the subscripts indicate the two- 
qubit subsystems. It is not difficult to show that the 
summations over all projectors \ej){ej\'^ (g) \ej){ej\ and 
\fk){fkV' ® \fk){fk\ in Eq. ^ factorize into products of 
n summations over two-qubit subsystems consisting of a 
single input and output qubit. The summations over the 
two-qubit subspaces can be performed with the use of the 
following identities 

|00)(00| -f |1I)(I1| = $+ -f 
I+ +)(+ + !+ !--)(--I = $+ + d/+, (B2) 


where = ^(|01)±|10)) are the other two Bell states 
and we used the notation $+ = |$+)(<I>+|. The identities 
(B2) allow us to rewrite the operator R as 

WRW^ =($+)®” _ ($+ + _ ($+ -(- 

-1- ($+-p $- -p -p 

(B3) 


where we used the four Bell states to express the identity 
on the two qubit Hilbert space as I = <I>+-|-<I>“-|-4'+-|-'I'”. 
Since W is unitary, operator WRW^ has the same eigen¬ 
values as R. Moreover, WRW^ is diagonal in the basis 
formed by tensor products of Bell states, hence the eigen¬ 
values can be directly determined from the expression 
(B3). All the projectors c onta ined in the first three terms 
of the right hand side of (B3) determine the 2"+^ — 1 di¬ 
mensional zero eigenvalue subspace and it is easy to see 
that on the remaining subspace the eigenvalue is one. 
Thus, all the eigenvalues of WRW"^ are nonnegative, 
which proves that i? is a positive semidefinite operator. 
This proves that the fidelity bound (13) holds also for 
the case, when the bases {|ej)}j=i, {\fk)}k=i are com¬ 
putational basis and its Hadamard transform on every 
qubit. 
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